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The puzzle of recently observed insulating phase of graphene at filling factor u = Q'm high magnetic 
field quantum Hlall (QH) experiments is investigated. We show that the magnetic field driven Peierls- 
type lattice distortion (due to the Landau level degeneracy) and random bond fluctuations compete 
with each other, resulting in a transition from a QH-metal state at relative low field to a QH-insulator 
state at high enough field at = 0. The critical field that separates QH-metal from QH-insulator 
depends on the bond fiuctuation. The picture explains well why the field required for observing the 
insulating phase is lower for a cleaner sample. 
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Introduction. — The intriguing quantum Hall effect 
(QHE) in graphene has attracted a lot of attentions in 
recent years 0, The quantum Hall (QH) plateaus 
was initially found to be [3] axy = (4e^//i)(n -I- 1/2) = 
ve^/h, n — 0, ±1,±2,.... The result can be understood 
in the same manner as the QHE in usual two-dimensional 
electron gases with graphene's special properties [1, 01: 
The factor of 4 comes from spin and valley (so-called 
K and K' points) degeneracy, and 1/2, in the terminol- 
ogy of edge states, is from the interesting splitting of 
n = Landau level (LL) at sample edges due to the 
particle- hole symmetry This QHE rule was soon 

violated in stronger fields, and anomalous plateaus of 
1/ = 0, ±1 [HI due to the breaking of spin and valley degen- 
eracy were discovered although the detailed symmetry- 
breaking mechanisms are still under debate The re- 
cent surprise comes from the discovery of insulating phase 
at charge neutral point (CNP) in a high magnetic field 
B by Checkelsky et. al vA because it disagrees with the 
early appealing theory [J| that explained well all exist- 
ing experiments then. This discovery was confirmed by 
the temperature-dependence study of magneto-transport 
near CNP (zero-energy) by Zhang et. al [8] . Experiments 
reveal following universal features. 1) Depending on 
the sample quality and applied field, graphene at CNP in 
QH regime can be either a QH-metal or a QH-insulator 
(see below). 2) A magnetic field can drive a graphene at 
CNP from the QH-metal state at relative low field into 
the QH-insulator state at high enough field. The critical 
field required for such a transition is lower for a higher 
quality sample. 3) The nature of the transition seem to 
be Kosterlitz-Thouless type (KT-type). 

Two possible n = states have been proposed and ex- 
amined. One is the QH-metal 0] that occurs when spin 
split is larger than the valley split as illustrated in Fig. 
1(a). The n — Q LL splits first into spin- up (blue curves) 
LL and spin-down (red curves) LL. The further splits 



of valley degeneracy near the sample edges create two 
counter-propagating charge and spin carrying edge chan- 
nels, which results in a residual conductance of 2e} jh. 
This is a topological state protected by the electron-hole 
symmetry The other possibility is the QH-insulator 
when the valley degeneracy lifts first with a smaller spin 
split primarily due to the Zeeman effect. As shown in Fig. 
1(b), the K—LL and K'—LL move symmetrically to the 
opposite direction of zero-energy (due to electron-hole 
symmetry). As a result, the topological edge channels 
are absent at zero energy, and a true insulator appears. 
The QH-insulator was previously excluded from a pos- 
sible n = state mainly due to the early experimental 
observation of saturated longitudinal resistance at low 
temperature |4|. However, the recent observation of in- 
sulating phase near CNP suggests that QH-insulator is 
also a possible zero-energy phase of graphene 0, Q ■ 

The mechanism of the breaking of fourfold degeneracy 
of n = LL has been discussed extensively in litera- 
ture Whether graphene at zero-energy in the 
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FIG. 1: (Color online) Schematic illustration of position- 
dependent energy of n = LL for QH-Metal (a) and QH- 
Insulator (b), respectively. The position (in unit of mag- 
netic length Ic = yJh/eB) is measured from the sample edge 
(j/o = 0) and eo represents energy gap between n = and 
n = 1 LLs. The blue and red lines represent the spin-up and 
spin-down states, respectively. 
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QH regime is described by QH-metal or QH-insulator 
has also been studied by Yang and Han [l5|. They sug- 
gested that a graphene will be a QH-metal if spin splits 
first while the graphene is a QH-insulator if valley splits 
first. However, the study did not address why the same 
sample can change from a QH-metal to a QH-insulator. 
Assume that there is a field-induced metal-to-insulator 
transition, Nomura, Ryu, and Lee [3| tried to identify 
possible order-parameter and to argue the transition be- 
ing the experimentally observed KT-type. In this letter, 
we would like to provide a plausible path that a graphene 
at CNP can change from a QH-metal state to a QH- 
insulator state as the magnetic field increases. The hope 
is not only to understand why both QH-metal and QH- 
insulator phases are possible in the same sample, but also 
to provide a guide for experimental manipulation of the 
transition between the two phases. We will show that 
bond fluctuations due to intrinsic ripples in a graphene 
can resist the K-K' split in low magnetic field, resulting 
in the QH-metal state at CNP. In the opposite scenario, 
one should expect a QH-insulator state. 

Picture. — It shall be useful to first present our pic- 
ture of why a transition from the QH-metal to the QH- 
insulator phases should be expected in graphene in a 
strong magnetic field. Firstly, LLs of graphene in a 
strong magnetic field are highly degenerated. Accord- 
ing to Peierls instability for a solid or Jahn- Teller effect 
for a molecule [l7| , the gra phe ne can lower its energy 
through a lattice distortion [l8[. The degree of distor- 
tion should be proportional to the magnetic field because 
the LL degeneracy is proportional to the field 13 , [3] . 
The distortion destroys K-K' (or sublattices A and B 
for n = LL) degeneracy. According to various estima- 
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the K-K' split is linear in B-field for 
an ordered distortion, similar to the usual Zeeman effect 
but bigger. Although there may be other mechanism of 
breaking vallev degeneracy such as electron-electron in- 
teraction d, [3, , but the Peierls instability should be 
robust and universal for highly degenerated systems like 
graphene that can also be regarded as a large molecule. 
Secondly, it is inevitable to get ride of ripples in graphene. 
The ripples create the intrinsic bond disorders, and the 
intrinsic bond fluctuations tend to suppress the K — K' 
split in n = LL because the bond fluctuation shall 
restore the inversion symmetry of A- and B-sublattices 
IS El. As a result, the lattice distortion and the in- 
trinsic bond fluctuations compete with each other. At 
relative low field, the fluctuations overtake the distor- 
tion effect, and there is no K-K' valley splitting. The 
graphene at the Dirac point is a QH-metal. However, 
at high enough field, the Peierls-type lattice distortion 
dominates, and K-K' splitting occurs. The Dirac point 
becomes a QH-insulator. 

Model and method. — Low energy excita- 

tions of graphene are from tt— electrons that 
can be modeled by a tight-binding Hamilto- 



6t,, 



\ia >< ja\ + h.c. 



tij is the 



<ij>,cr 

hopping energy between two nearest-neighbor sites whose 
value for a pure graphene is i = —2.7eV. ct = ± is spin 
label. is the on-site energy that includes the Zeeman 
contribution ±gfiBB /2 for spin-up (+) and spin-down 
(— ), leading to a spin-gap As = g^J^sB — 1.3 x B[K/T]. 
One can also add a random on-site energy to mimic 
the extrinsic disorders. 5tij is for the possible intrinsic 
random hopping energy described later. The magnetic 
field is introduced by Peierls' substitution in hopping 
parameter 4>ij — 2TTe/h J:' A ■ dl The Peierls lattice 
distortion is either in-plane like Kekule bond order 
which alternates the short and long 
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bonds like in the classical benzene molecule (see inset of 
Fig . 3(a)), or out-plane like charge density wave (CDW) 
[l8| , which breaks the inversion symmetry of sublattices 
(see inset of Fig. 3(d)). The distortion is chosen in 
such as way that valley gap is = 2.0 x B\K/T] 
EH [Hi (for in-plane like Kekule bond order wave) 
or A„ = 4.2 X B[K/T] [3 (for out-plane CDW), both 
of which are bigger than A^. 

It is known that ripples (long-ranged corrugation) are 
the intrinsic fiuctuations of graphene (2^ . It is observed 
that the graphene ripples have a height variation of about 
0.5 ~ l.Onm, which occurs on a lateral scale of 8 ~ lOnm. 
The graphene with random ripples can be simulated by 
superposition of a series of plane waves [2^. The out- 
of-plane displacement is h{r) = C^C^. sin(qj • r -I- 5i), 

i 

where r = {x,y) is the position of in-plane atoms and 



Cq . , q.j and 5i are random numbers 25| . The lowest pos- 



sible Qj determines lateral scale of ripples, and the height 
fluctuation is controlled by C value. Fig. 2 is the contour 
plot of one corrugated graphene used in our calculations. 
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FIG. 2: (Color online) Contour plot of a randomly generated 
landscape of the ripples used in our calculations. The average 
lateral scale of ripple is about Xrippie = 8nm. The height 
fluctuation is approximately in the range of [-1.0,1.0]nm. 
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The random shapes and locations of ripples lead to a ran- 
dom change in bond length that in turn causes the hop- 
ping energy fluctuations [2|: Stij — taAa/a, where Aa = 
+ (/i(r) - /i(r'))2-a and a = dlogt/dloga = 2^. 

As it was explained in our early discussions, whether a 
graphene is in QH-metal or QH-insulator states depends 
on how the valley and spin split in the bulk. Thus one 
needs only to study bulk states in order to determine two 
topologically different states. To see how n = LL split 
into spin bands and valley bands in the lattice distorted 
graphene with intrinsic random hopping, we need to ob- 
tain an accurate density of states (DOS) of the model 
numerically. The Lanczos recursive method [2^ [13] on a 
large lattice (with about million lattice sites) is employed. 
The averaged DOS p{E) ^ < -0| e-h+iti \^ > 

computed. In the approach, a small artificial cut-off en- 
ergy (77 = O.lmeV), which results in a small LL width in 
clean graphene, is introduced to simulate an infinitesimal 
imaginary energy p7| . A periodic boundary condition is 
used in order to minimize the boundary effects since our 
concern is the bulk n — LL splits. 

Results and discussions. — Fig. 3(a-c) are the DOS of 
n = LL for in-plane distorted graphene of the Kekule 
bond order wave (shown in the inset of Fig. 3(a)) in 
the presence of ripples shown in Fig. 2. At a relative 
low magnetic field B = 9.4r (Fig. 3(a)), the distortion- 
induced valley gap cannot compete with the hopping fiuc- 
tuation so that bond order wave is statistically destroyed. 
As a result, K-K' degeneracy is statistically restored, 
and n = LL splits into spin-up (blue) and spin-down 
(red) bands due to the Zeeman energy as shown in Fig. 
3(a). This is the QH-metal state discussed early because 
the spin splitting dominates the bulk gap of n = LL. As 
the magnetic field is above the critical field Be (between 
9.4T and 13. 2T in our case as indicated in Fig. 3(a-b)), 
the distortion-induced valley gap appears and n = LL 
splits into four bands. As shown in Fig. 3(b), two peaks 
on the left (right) sides of zero energy are for K (K') val- 
leys and the corresponding edge states do no cross each 
other. The graphene is in the QH-insulator state at CNF. 
Further increase of the field, the valley splits become even 
bigger as shown in Fig. 3(c) for 21. 9T. To shown that 
above result is very robust against different type of dis- 
tortion, we made the same type of calculations for the 
CDW lattice distortion (shown in the inset of Fig. 3(d)) 
[Tsf . The similar results as those for the Kekule bond 
wave have been found. As shown in Fig. 3(d-f). The 
valley degeneracy is preserved at 6.6T (d) and the zero- 
energy state is a QH-metal. The valley degeneracy is 
broken at 11. OT (e), resulting in a QH-insulator at CNF. 
Further increase of field to 21.9T (f) in this case, the two 
K [K') bands are clearly located in the same side of the 
zero-energy point. 

Fig. 3 shows clearly that the intrinsic bond random- 
ness due to ripples can suppress lattice distortion-induced 
valley splitting no matter whether the distortion is in- 
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FIG. 3: (Color online) Density of state of n = LL of lat- 
tice distorted graphene in the presence of intrinsic random 
hopping for various magnetic field, (a-c) are for in-plane dis- 
tortion of Kekule bond wave in B =9.4T (a); 13. 2T (b); and 
21. 9T (c). (d-f) are for out-plane CDW with B =6.6T (d); 
IIT (e); and 21. 9T (f). The blue (red) line represents the 
spin-up (spin-down) state. 



plane (Kekule bond wave) or out-plane (CDW), demon- 
strating the robustness and universality of the picture 
for grapheme. As shown from our calculations, the val- 
ley degeneracy is destroyed and valley gap appear when 
the bond disorder is not strong enough to restore inver- 
sion symmetry of A-B sublattices statistically. A cleaner 
sample should have a weak bond disorder, and the cor- 
responding valley degeneracy can be broken by a smaller 
lattice distortion. In another word, the critical field Be 
that induces the transition from the QH-metal state to 
the QH-insulator state depends on the degree of bond 
randomness. The cleaner a sample is, the smaller the 
critical field required to drive the graphene from the QH- 
metal into the QH-insulator is needed. In our model the 
bond randomness is measured by the height-fluctuation 
hmax of ripples. Fig. 4 is /imax-dependence of the criti- 
cal field Be for both Kekule type and CDW-type lattice 
distortion. It is found that the value of Be of CDW type 
is lower than that of Kekule type because the valley gap 
in CDW distortion is bigger. The inset of Fig. 4 is the 
phase diagram for the QH-insulator and the QH-metal. 
This is exactly what was observed in experiments 0, 0] ■ 
This also explains that failure of the discovery of the 
zero-energy insulating phase at high field in experiments 
before Checkelsky et. al is probably due to the quality 
of the samples. 

We would like to make a few remarks before ending this 
paper. 1) Fixed the applied magnetic field, the random 
fiuctuation drives a graphene from a zero-energy insula- 
tor to a zero-energy metal. This is opposite to the role of 
randomness usually played in a metal-to-insulator tran- 
sition. 2) There are two types of disorders in graphene, 
the intrinsic bond disorder due to the spontaneous rip- 
ple formation and extrinsic disorder from the impurity 
states in substrate (SsJ and/or ad-atoms The ex- 
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FIG. 4: (Color online) Relationship between critical field Be 
and ripple effects for Kefcule bond order (magenta) and CDW 
(olive), respectively, hmax represents the maximum height 
fluctuation of spatial ripples of samples. Inset: Phase diagram 
of the QH-metal state and QH-insulator state. 



trinsic disorder can cause the on-site-energy random- 
ness in the tight-binding model. Our numerical studies 
show that the on-site disorders will not suppress the K- 
K' split, but merely broaden K— and K'— subbands 



12l [29|. The Landau subband broadening can cause 
inter-Landau-subband mixing that may change the lo- 
calization property [s^ of a disordered system. 3) Our 
QH-metal is in fact a topological insulator because it is 
a bulk insulator with conducting edge channels Q while 
the QH-insulator is a conventional insulator. The present 
work is on why the magnetic field can drive such a topo- 
logical transition rather than the nature of the transi- 
tion [l6|. 4) The Coulomb interaction is neglected here 
because we do not expect the change of physics by the 
interaction. 

We present a theory of graphene for the field-driven 
topological transition from the QH-metal state to the 
QH-insulator state at = 0. The transition results from 
the competition between Peierls-type lattice distortion 
and random bond fluctuation. Our theory provides not 
only a clear explanation about the existence of both QH- 
metal state at low field and QH-insulator state at high 
field near the CNP, but also why the critical field for a 
cleaner sample is lower. 
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